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Abstract 

The purpose of this paper is to develop the negative order MKdV hierarchy 
and to present a new related integrable Neumann-like Hamiltonian flow from 
the view point of inverse recursion operator and constraint method. The 
whole MKdV hierarchy both positive and negative is generated by the kernel 
elements of Lenard's operators pair and recursion operator. Through solving 
a key operator equation, the whole MKdV hierarchy is shown to have the 
Lax representation. In particular, some new integrable equation together 
with the Liouville equations, the sine-Gordon equation, and the sinh-Gordon 
equation are derived from the negative order MKdV hierarchy. It is very 
interesting that the restricted flow, corresponding to the negative order MKdV 
hierarchy, is just a new kind of Neumann-like system. This new Neumann- 
like system is obtained through restricting the MKdV spectral problem onto 
a symplectic submanifold and is proven to be completely integrable under 
the Dirac-Poisson bracket, which we define on the symplectic submanifold. 
Finally, with the help of the constraint between the Neumann-like system and 
the negative order MKdV hierarchy, all equations in the hierarchy are proven 
to have the parametric representations of solutions. In particular, we obtain 
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the parametric solutions of the sine-Gordon equation and the sinh-Gordon 
equation. 
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1 Introduction 

It is well-known that the nonlinear evolution equations (NLEEs) solved by the fa- 
mous inverse scattering transformation (1ST) can be understood as compatible con- 
ditions of some linear equations |1|, namely. Lax representation. In the past two 
decades the Lax representation has played a very important role in the discussion 
of NLEEs. In particular, the Lax representation has been used successively in the 
bi-Hamiltonian structure of finite-dimensional dynamical systems [^, |^ , in the non- 
linearization theory of soliton system to produce new completely integrable systems 



in the Liouville sense |Tl], |12] , in the tri-Hamiltonian formulation of nonlinear equa- 
tions |T^, and in the finite-dimensional restricted flows of the underlying infinite 
systems 0. Recently, the finding of peaked solitons produced a breakthrough in 
the study of nonlinear partial differential equations [Q, where Camassa and Holm 
showed their equation is completely integrable by the 1ST method. This fact allows 
one to discuss the nonlinear dynamics of soliton solutions and billiard solution 0] via 
their linear spectral content (i.e. Lax representation). Thus, for a given hierarchy 
of NLEEs, to find the Lax representation is of great importance. 

The modified Korteweg-de Vries (MKdV) equation together with the MKdV hi- 
erarchy has wide applications in physics and other nonlinear sciences. It possesses 



the Lax pair the periodic soliton solution ||3y, the bi-Hamiltonian structure 
[0 and other soliton properties such as Darboux transformation, Backlund transfor- 
mation and the Miura transformation between it and the KdV equation pOl . About 



the study of the MKdV hierarchy, i.e. usual higher-order MKdV equations, there 



have been many discussions in the literature. In Ref. |15[, the sine-Gordon equation 
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and the intergrated MKdV equations were shown to conserve the same infinite set of 
charges which were determined by a recursion relation. Afterward, Verosky p9 



in- 



troduced the negative powers of Olver recursion operator and presented the relations 
between the sine-Gordon/sinh-Gordon equation and the potential MKdV equations 
(non-local fiows). In 1993, Andree and Shmakova discussed the supersymmetry 
structure of the sine-Gordon equation, which can be embraced in the MKdV and 
KdV hierarchies. All of those interesting facts happened between the sine-Gordon 
equation and the MKdV equation. The fact that the sine-Gordon equation is a 



negative fiow of MKdV can be seen in Ref. 



In recent years, the time-discrete version of integrable systems have already 
arisen a lot of attractive attentions. This idea was applied to constructing a new Lax 
pair from the old Lax pair ||22| . In Ref. ||2^ , we proposed an approach to generate the 
positive and the negative order integrable hierarchies from a given spectral problem. 
In fact, we had this idea in Ref. [^. But that depends on the existence of inverse 
recursion operator. This means, for every concrete spectral problem, we need to 
determine the inverse of recursion operator. 

In the present paper, we will give the inverse of recursion operator of the MKdV 
hierarchy in an explicit form. With the help of the recursion operator and its inverse, 
we present the positive order and the negative order MKdV hierarchies of NLEEs. 
The whole MKdV hierarchy is proven to have the Lax pair through employing a key 
operator equation, and is therefore a completely integrable hierarchy. Particularly, 
some new integrable equation together with the Liouville equations, the sine-Gordon 
equation, and the sinh-Gordon equation are derived from the negative order MKdV 
hierarchy. Furthermore, the constraint between the MKdV spectral problem and the 
negative order MKdV hierarchy derives a restricted Hamiltonian fiow on a symplectic 
submanifold, which is a new kind of Neumann-like system. Then, we define the 
Dirac-Poisson bracket on the symplectic submanifold, which is a very useful tool to 
deal with the finite dimensional integrable system on some submanifols. Under the 
Dirac-Poisson bracket, the new Neumann- like system has the Hamiltonian canonical 
form and has furthermore independent and involutive functions, which guarantees its 
complete integrability. Finally, each equation in the negative order MKdV hierarchy 
is proven to have the parametric representation of solution, which is given by the 
involutive solution of Hamiltonian phase fiows (i.e. x-fiow and t^-fiow, m < 0, m G 
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Z). Particularly, we obtain the parametric solutions of the sine-Gordon equation 
and the sinh-Gordon equation. 

The whole paper is organized as follows. Next section gives the pair of Lenard's 
operators, the recursion operator and their inverses, and introduces a key operator 
equation which is available for both the positive order and the negative order MKdV 
hierarchies. This operator equation is different from the one usually considered in 



literatures [|T0|, ^Tj. In section 3, we will see that the positive order MKdV 
hierarchy, arising from the kernel of one of Lenard's operators, is nothing but the 
well-known MKdV hierarchy. In section 4, we present the negative order MKdV 
hierarchy by using the kernel element of the other one of Lenard's operators. All 
equations in the hierarchy have the Lax pairs. In sections 5 and 6, we provide a 
new kind of integrable Neumann-like system on a symplectic submanifold, and give 
the parametric solutions of the negative order MKdV hierarchy, respectively. In the 
last section, we give some conclusions. 

For convenience, we make the following conventions: 

Qx-ni f fmx ; 0,1,2,..., 

/(-) = <! / fdx, m = -l,-2,.... 



ft = %, fmxt = 15^^ ("^ = 0, 1, 2, ...), 9 = ^, (9 Ms the inverse of d, i.e. dd ^ 
d~^d = 1, 9™"/ means the operator d"^f acts on some function g, i.e. 

r £,ifg) = ifgU., m = 0,1,2..., 

d'^f-g = d'^ifg) = \ J ...J fgdx, m = -1, -2, .... 



In the following the function u stands for potential, the imaginary unit i is 
satisfying P = —1, and A is assumed to be a spectral parameter, and the domain of 
the spatial variable x is which becomes equal to {—oo, +oo) or (0, T), while the 
domain of the time variable tm is the positive time axis = {tm\ tm G M, tm > 
0, m = 0, ±1, ±2, ...}. In the case Q = {—oo, +oo), the decaying condition at 
infinity and in the case fl = (0, T), the periodicity condition for the potential 
function is imposed. 
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(M^^,(ip A dq) stands for the standard symplectic structure in Euclid space 
M2iv = I (^p^ p = (p^^ . . . ^p^)^ g = (g^^ . . . ^ g^)}^ p^.^ g^. (j = 1, . . . , AT) are pairs 

of canonical coordinates, (■, ■) is the standard inner product in M^; in (M^^, dpAdq), 
the Poisson bracket of two Hamiltonian functions F, H is defined by 



j^^\dqjdpj dpjdqj 



N ' d£ dH\ _ / d£ dH_' ^ ^^^^ 

dq^ dp / \ dp^ dq 



Xi, Xn are distinct spectral parameters, and A = diag{Xi, Xn)- Denote all 
infinitely times differentiable functions on real field M and all integers by C°° 
and by Z, respectively. 



2 Inverse recursion operator and operator equa- 
tion 

Let us consider the following spectral problem 

which is quite a special case of the well-known Zakharov-Shabat-AKNS spectral 



problem ^| 



^.= 1 "'^ 1^ (2.2) 
\ V iX 

with V = u. 

Eq. (|2.1|) is equivalent to 

L.^^ I ''1 1 ■^ = A^ (2.3) 

lu —to 

and its spectral gradients VA = |^ = — -0^ satisfies the Lenard eigenvalue problem 

K ■VX = X^J -VX (2.4) 
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with the Lenard's operators pair 



K = --d^ + dud-^ud, J = d, 
4 



which yields the recursion operator 



1 



£ = J-^K = — + ud~\d. 
4 



(2.5) 



(2.6) 



Apparently, the Gateaux derivative operator L^,((^) of the spectral operator L 
given by Eq. (O) in the direction ^ G C°°(R) is 



L{u + ee) 



e=0 



-i^ 





(2.7) 



which is obviously an injective homomorphism. 

Through guesswork and calculations, we can obtain the inverse operators of 
L, J, K and £: 



d~^uA -A 



d 



-1 



(2.8) 

(2.9) 
(2.10) 
(2.11) 



For any given C°°-function G, we construct the following operator equation with 
respect to = V{G) 



[V,L] = L,{K ■ G) - L,{J ■ G)L^ 



(2.12) 



Remark 1 This equation contains a special term instead of the term L usually 



considered in literatures / |7^ \2^, [F^/ , 
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Theorem 1 For the MKdV spectral operator ( ^.Sj ) and an arbitrarily given C° 
function G, the operator equation ^2.12^ ) has the following solution 

T/ _ — I (^^^)*' ^^d—^uGx 2^xd — jGxx i /oiq^ 



Proof: Directly substituting Eqs. (|27[3|) , and into Eq. (gl2D , we 

can complete the proof of this theorem. 



3 Positive order hierarchy of Eq. ( |27ll ), i.e. usual 
MKdV hierarchy 

Let us now give the positive order MKdV hierarchy through considering the kernel 
element of the Lenard's operator J. 

Gq = a & KerJ and the recursion operator (|2.6| ) yield the positive order hierar- 
chy of Eq. (O) 



ut,^ = JC"" -a, m = 0,1,2,..., (3.1) 
which has the following representative equations 

Ut^ = aux, trivial case, (3-2) 
= -^auxxx + ^au^Ux. (3.3) 

Here, a = a(t„) G C°°(]R) is an arbitrarily given fuction with respect to variables 
tn (n > 0, n G Z), but independent of x. Apparently, with a = 4 Eq. (|3.3| ) becomes 
the well-known MKdV equation 

- Qv?Ux + Uxxx = 0. (3.4) 

Therefore, Eq. ( p.l|) coincides with the well-known MKdV hierarchy, which corre- 
sponds to the isospectral case: Xt^ = 0. 
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By Eq. (|2.13| ), the whole hierarchy has the standard Lax representation 

= [iy^,L], (3.5) 



771—1 



= J2^{Gj)L'^"'~^~'\ (3.6) 

j=0 

where V{Gj) is given by Eq. ( p. 131) with G = Gj = D - a, j > 0, j G Z. Therefore 
we obtain the following theorem. 

Theorem 2 The positive order hierarchy (i.e. the MKdV hierarchy) of the 

spectral problem ( ^.dj ) possesses the Lax pair 

/-■X \ ^^-^^ 
Ar.=a{ 1 A2(-^V + Eri'^.^'^'""^"'V' m = 0,l,2,..., 



where 



(3.8) 

Gj = a - a, j> 0, j e Z. (3.9) 



In particular, the MKdV equation ( p.3| ) has the Lax pair 




- ^ +1^.3 ^(A^ + iAn^^ 



Remark 2 For the spectral problem ^2. if we use the usual method pS^ , ^q] , i.e. 
the method of finite power expansion with respect to spectral parameter X, then no 
isospectral evolution equations of Eq. (\2. H) can be obtained. 
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However, we here present the MKdV hierarchy (\3. 1\ ) purely by the Lenard's op- 
erators pair satisfying Eq. ( ^.4 )- Due to containing the spectral gradient VA in Eq. 



(\2.^ ), this procedure of generating evolution equations from a given spectral problem 



is called the spectral gradient method. 

In this method, how to determine a pair of Lenard's operators associated with the 
given spectral problem mainly depends on the concrete forms of spectral problems 
and spectral gradients, and some computational techniques. From this method, 
we can furthermore derive the negative order MKdV hierarchy, which is displayed 
below. 

4 The negative order MKdV hierarchy and Lax 
representation 

Let us now give the negative order MKdV hierarchy through considering the kernel 
element of Lenard's operator K. The kernel of operator K has the following three 
seed functions: 

G\ = fe--'^^0'"^ (4.1) 



G\ = (e-'-^')^"\ (4.2) 
Gt, = (e-'-'-''y~'Ue'-'-'')^~'\ (4.3) 



whose all possible linear combinations form the whole kernel of K. Let G_i G 
Ker K, then 

3 

G_i = 5^a,G'^ (4.4) 

k=l 

where = afc(t„), k = 1,2, 3, are three arbitrarily given C°°-functions with respect 
to variables t^ {n < 0,n G Z), but independent of x. Therefore, G^i directly 
generates the isospectral (A^^ = 0) negative order hierarchy of nonlinear evolution 
equations for the spectral problem ( |2.3|) 

ut^ = JC^^^ ■ m < 0, m G Z, (4.5) 
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which is called the negative order MKdV hierarchy of Eq. ( |2.3| ). By 

Theorem |l], the hierarchy ( [4. 51) has the standard Lax representation 

= [iy^,L] (4.6) 
-1 

= -Y,V{Gj)L'^"'-'-'\ m = -l,-2,..., (4.7) 



I.e. 



with 



= ' u Ta ) ^' (4.8) 

= EjL 1^, A2(™~^-~i)^, m = -1, -2, 



where 



C ■ Gj — —-Gj^xx + u {'^Gj^xY '* 



In Eq. (13), 1/ {Gj) and L'^ are given by Eq. (^) with G = = ■ G_i 
and by Eq. (|2.^ ), respectively. Thus, all nonlinear equations in the hierarchy (f4.5|) 
are integrable. 

Let us now give some special reductions of Eq. ( ^.5| ). 

• In the cases of 02 = cts = 0; ai = = 0; ci = 02 = 0, Eq. ([4.5|) separately 
has the following representative equations 

(4.10) 
(4.11) 

2„(-. ^,2«(-)y-^) ^ ^2«(-) (^e-2«<-'y-^)^ , (4.12) 

which can be via the transformation u = Vx respectively changed to 

■J^x.i.i = aie~'^^ , Liouville equation, (4-13) 
Vx,t.i = 0.26'^'", Liouville equation, (4-14) 
Vt-i = 0-3 {g^^Y ""^^ (e"^^)'' ^\ a new integrable equation. (4.15) 
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They are also equivalent to the following differential equations: 



Ux,t.i + 'iuut^^ = 0, 
Ux,t^i - '2uut_^ = 0, 

Uxx,t^i - U~^U^,Ux,t_i + U'^Ux - 4:u'^Ut_^ = 0. 



(4.16) 
(4.17) 
(4.18) 



Apparently, Eqs. ( [4.10| ), ( [4.11| ), and ( [4.12| ) possess the following standard Lax 
pairs 



i^x 



-iX u 
u iX 



i^x 



-iX u 



u iX 
-iX u 



1 -1 



u iX 



1 

1 



A-^V^ + [/_iA-V, 



respectively, where 



(4.19) 



(4.20) 



(4.21) 



(4.22) 



Eqs. ( [4.16| ) and ( |4.17| ) are two Liouville equations, and easy to see that they 
can be directly integrated as u'^ ±Ux = f{x), f{x) G C°°(]R), which are two typical 
Riccati equations. They can be solved by some methods in the theory of ordinary 
differential equations. But, Eq. ( |4.15| ) or ( [4.1^ ) is a new integrable equation. 



11 



In the case of ai 
reads 



\, 0'2 = :|, and = 0, the first equation of Eq. ([4 .51) 



We make a simple transformation 



U = -IVx 

2 



(4.23) 



(4.24) 



then Eq. ( ^.23| ) is exactly changed to the well-known sine-Gordon equation 

^x,t_i = sin v. (4.25) 



According to Eq. 
pair: 



the sine-Gordon equation ( |4.25|) possesses the following Lax 



-iX h' 



2 '^'^x 



iX 



_1_ 

4A 



I COS V sm V 
— sin V —i cos v 



(4.26) 



which has a slight difference from the usual one given in Ref. [0, [17] 
For Eq. ( [4.23| ), if we make the transformation u 



^Vx, then it becomes 



Vx,t^i = sinh V 



(4.27) 



which is nothing but the well-known sinh-Gordon equation. By Eqs. ( [4.8|) and 



the sinh-Gordon equation (|4.27|) has the following Lax pair 



-iX \v.j, 

:Vr iX 



1^, 



cosh V 
sinh V 



sinh V 
coshf 



(4.28) 



which is also slightly different from the usual one ||3^, [17 



Remark 3 In Ref. j^, Alber, Camassa, Holm and Marsden took evolution equa- 
tions of auxiliary linear system polynomial in X^^ and gave the zero- curvature rep- 
resentation for the Dym type hierarchy and the Camassa-Holm equation. Here, our 
starting point is the spectral problem and the spectral gradient instead of the auxiliary 
linear problem, then via a key operator equation we obtain the Lax pairs. 
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5 A new integrable Neumann-like system 



In Ref . , we presented an integrable Neumann-like system closely associated with 
the positive order MKdV hierarchy (|3.1| ), showed that the Neumann-like system was 
sent by a gauge transformation to an integrable Bargmann system, and obtained 
the parametric solutions for the positive order MKdV hierarchy. Now, we study the 
negative case. 

Let Ai, Aiv be different spectral parameters of Eq. (|2.1|), {qj,Pj)^ the spec- 
tral function corresponding to Xj, and P = —1. Then Eq. becomes 

qj,x = -iXjQj + upj, 

Pj,x = UQj + iXjPj, 

i.e. Eq. (O) yields the following form: 

= -iAq + up, 
Px = uq + iAp, 

where p= {pi, ...,pn)^ , q= (qi, q^)^ , and A = diag{\i, Xn)- 

Assume {Aq,p) 7^ 0, then let us restrict Eq. ( pTip onto the following symplectic 
submanifold M in M^^ 



{q,p) e 



t,2N 



G^^{{Aq,q) + {Ap,p)) = 0}. (5.2) 



Then, on M we obtain the following constraint 

^ (AV,>-(AV.ri 
2{Aq,p) 

Thus, under Eq. ( |5.3| ) the MKdV spectral problem (|2.1|) is nonlinearized as the 
following nonlinear system 

(A^g,g)-(A^p,p)^ 



g,. = -iAq + i^-^ix^^p^P, 
= ^- — — -Q + ^^P^ 



2(Ag,p) 

which is called a restricted MKdV flow of the spectral problem ( |2.1|) on 
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Remark 4 Because dF, dG are independent everywhere on M and their deter- 
minant det({F, G}) = {Aq,p) ^ 0, M. is therefore a symplectic submanifold in 
/|7?|/. Apparently, M is not the usual tangent bundle, i.e. M 7^ TS^^^ = 
I {q,p) e M^^l F = {q,q) - 1 = 0,G = {q,p) = o|, thus Eq. does not coin- 

cides with the usual Neumann constraint / [7^ \TBJ on TS^^^. If we strongly impose 
Eq. I\5.1\ ) on the usual tangent bundle TS^~^ = | {q,p) G M^^| F = g) — 1 = 0, 

G = {q,p) = 0}, then we have no constraints except for u = which is of course 
meaningless. So, Eq. I \5.^ ) has no link to the standard Neumann system and is 
therefore a new kind of Neumann-like system. 

In order to prove the integrability of the restricted flow ( |5.4| ) on M, we introduce 
the Dirac bracket 

{/, 9}d = {/, g} + ({/, F}{G, g} - {/, G}{F, g}) (5.5) 

which is easily proven to be bilinear, skew-symmetric and satisfy the Jacobi identity. 
Let us now consider a very simple Hamiltonian function 

H = -i{Aq,p) (5.6) 

together with independent functions 

= l{{A'-p,p)-{A'-q,q)) 

+\ E - {A'^'^%,P)) {{A'^-~'-%,p) - (A^(-^-"^)g,g» 

j=m 



1^ 



4 

j=m 



(A^^'+ig, q) + (A^i+V, p) 2 (A2(™-^')- V, q) 

2 (A^J'+ig, p) (A2("^-J)-ig, q) + (A2(™-J)-ip^ 

m = — 1, —2, .... 

(5.7) 



Lemma 1 The inner product ^§^1 symmetric with respect to m, n (m,n < 

0,m, n E Z). 
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Proof: Making the derivatives of Fm with respect to q, p and directly substituting 
them into have a lengthy calculation and then know that this inner 

product is sum of some symmetric terms with respect to m, n (m, n = —1, —2, ...). 

Proposition 1 

{F,„, F„} = {H, E^} = 0, m, n = -1, -2, .... (5.8) 

Proof: Lemma |I| directly yields {Fm, Fn} = 0. As for the second equality, a 
straightforward computation completes its proof. 



Through some guesswork, we find that the restricted fiow (|5.4| ) can be expressed 
as the canonical Hamiltonian form in the following theorem. 



Theorem 3 Under the Dimc-Poisson bracket the restricted MKdV flow ( \5.^) 

coincides with: 

{H)d ■■ < r (5.9) 

[ Px = {P,H\d, 

where H is defined by Eq. ( |5.6[ j. 

A furthermore direct calculation leads to the following lemma. 
Lemma 2 

{Fm,Fn]D = {H,Fm}D = ^, m, 71 = -1,-2,.... (5.10) 
Because Fm are independent, we obtain the following theorem. 

Theorem 4 The restricted MKdV flow on the symplectic submanifold M is 

completely integrable. Moreover, all restricted flows {Fm) /? on M 

{Fm)D: I ^*'" = |^'^->^' (5.11) 

are integrable. 
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6 Parametric solution of the negative order 
MKdV hierarchy 

In the following, we will consider the relation between constraint and nonlinear equa- 
tions in the negative order MKdV hierarchy ( ^.51 ). Let us start from the following 
setting 

N 

^-1 = EV^^^' (6-1) 

i=i 

where G-i are defined by Eq. (|4.4| ), and VAj = Pj — q] is the functional gradient of 
the spectral problem ( p.lD corresponding to Xj. 

Let the recursion operator C act on both sides of Eq. ( [OD . Then, through 
making a choice of ■ = d'^ ■ = ^, we obtain 

{p,p) - {q,q) = (6.2) 



Doing the derivative on both sides of Eq. ( |6.2| ) with respect to x and substituting 
Eq. yield 

{Aq,q) + {Ap,p) = 0, (6.3) 



which together with Eq. ( |6.2| ) forms the symplectic submanifold M we need. Ap- 
parently, derivative for Eq. ( |6.3|) with respect to x leads to the constraint relation 

(U- 

Since the restricted Hamiltonian flows {H)ij and {F^jD ai'e completely inte- 
grable and their Poisson brackets {H, Fm}D = (m = —1, —2, ...), their phase flows 
gf[, g^p commute Thus, we can define their compatible solution as follows: 



q{x,trn) \ ^ . ( qix°,t'^ 
p{x,tm) / I p{x',t 



9h9fZ{ :)^o,7i I' rn = -l,-2,..., (6.4) 



where t'^ are the initial values of phase flows gfj, g^p^. 

Theorem 5 Let q{x,tm), p{x,tm) be a solution of the compatible Hamiltonian sys- 
tems {H)£) and {Fm)D on M. Then 

.{A'^q{x,tm),q{x,tm)) - (A^p(x, t„),p(2;, t„)) 
u = I — — — (6-5) 

2 {Aq{x,tm),p{x,tra)) 
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satisfies the negative order MKdV hierarchy 

ut^ = JC'^+'-G^,, m = -l,-2,.... (6.6) 



Proof: On one hand, the recursion operator C acts on Eq. (|6.1| ) and resuhs in the 
following 

j^m+l ^ J . A2>, p) - (A2™g, g)) 

= 2((A2>,p,)-(A2-g,g,)) 

= 22((A2"+V,p> + (A'"^+^?,g». (6.7) 



In this procedure, Eqs. (|2.4|) and (|5.4|) are used. 

On the other hand, we directly make the derivative of Eq. (|6.5|) with respect to 
tm- Then we obtain 

2(Ag,p) 

- {{A\,q) - {A'p,p)) {{Aq,PtJ + {Ap,qtJ)) (6.8) 
where q = q{x, tm), p = p{x, tm). But, 

d F^Yfl d ^TfL / \ 



therefore after substituting them into Eq. (|6.8|) and calculating it, we have 

ut^ = 2i ((A2™+ip,p) + (A2'"+ig,g» (6.10) 
which completes the proof. 

In the special case of m = —1, we have the following corollary. 

Corollary 1 Let q{x,t^i), p(x,t_i) be a solution of the compatible Hamiltonian 
systems {H)o and {F^i)£) on M. Then 

.(A2g(x,t_i),g(x,t_i)) - {A'^p{x,t-i),p{x,t_i)) 

U = % T-, -, ^ 7 TT (6.11) 

2(Ag(x,t_i),p(x,t_i)) ^ ^ 
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is a solution of the nonlinear evolution equation Ut_^ = X]fc=i '^'=(^-1)^-1 

+a, ( e''-'-' (e'-'-'Y''^ + e^"^"^' f e-^"^"^^) ^"'^^ , (6.12 



Here H is defined by Eq. ^5. dj ) and F_i is given by 

~ {A-\q + p),q + p) {A-\q -p),q-p). (6.13) 

In particular, the Liouville euqations v^^t-i = aie~'^'" , v^^t-i = 0,2^'^^ i the sine-Gordon 
equation v^^t-i = sinv, and the sink-Gordon equation v^^t-i = sinhf have the fol- 
lowing parametric solution 

(_i) . /■ (A2g(x,t_i),g(x,t_i)) - (A2p(a:,t_i),p(x,t_i)) 

V = '=11 — ; : ; dX , (6.14) 

J 2{Aq{x,t.i),p{x,t.,)) 

V = -2.n(-) = / (AM^^t_,),g(x t,,)) - (AMx t i),p(x,t_0) 

J {Aq{x,t.i),p{x,t^i)) ^ ' 

and 

_ o,,(-i) _ , [ (A'g(x,t-i),g(x,t_i)) - (AVx,t,i),p(x,t-i)) 

~ J {Aqix,t_,),pix,t.,)) ^'-^'^ 

respectively. 

By Theorem ^ the constraint (|5.3| ) is a solution of the negative order MKdV 
hierarchy ( |6.6| ). Thus, we also call the system {H)d (i.e. Eq. ( |5.4| )) a negative 
order restricted MKdV flov^^ of the spectral problem ( p.l|) on the sym- 
plectic submanifold M. All Hamiltonian systems (Fm)/^ (i.e. Eq. (|5.11| ) 
derived from {H)]j) are therefore called the negative order restricted flows 
on M. 



7 Conclusion and Comparison 

It is well-known that some traveling wave solutions or soliton solutions for the in- 
tegrable equations are available by the Inverse Scattering Transformation. Thus, a 
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natural question is: what is the relationship between the traveling wave solutions 
and the parametric solutions ( |6.5|) for the MKdV case? Because we need either the 
potential function u decaying at ±00 or satisfying the periodic condition in the pe- 
riod T with respect to the variable x (see the part of Introduction) when we do the 
spectral gradient calculations, we believe that both the traveling wave solutions and 
the periodic or quasi-periodic solutions should be in the formula (|6.5|) , namely, they 
share a common expression (|6.5|) . In a further procedure, we will consider giving an 
explicit expression of Eq. ( |6.5| ) for the cases of the potential u decaying at ±00 or 
satisfying the periodic condition in the period T. 

For the positive order constrained MKdV system by the constraint u = {p, p) — 
{q,q), we have dealt with it in Ref. where we knew this constraint is closely 
connected to the positive order (i.e. usual) MKdV hierarchy ( |3.1|) and in detail 
discussed the integrability of the constrained flow for the spectral problem (PTT|). 

A systematic approach to generate new integrable negative order hierarchies of 
NLEEs can be seen in Ref. pll . 
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